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Abstract 

We have applied a new noncompact, gauge- invariant, Monte Carlo 
method to simulate the U(l), SU(2), and SU(3) gauge theories on 
8 and 12 lattices. For U(l) the Creutz ratios of the Wilson loops 
agree with the exact results for (5 > 0.5 after a renormalization of the 
charge. The SU(2) and SU(3) Creutz ratios robustly display quark 
confinement at /3 = 0.5 and (3 = 2, respectively. At much weaker 
coupling, the SU(2) and SU(3) Creutz ratios agree with perturbation 
theory after a renormalization of the coupling constant. For SU{3) 
without quarks, our lattice QCD parameter is = 130 ± 18 MeV. 
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1 Introduction 



In compact lattice gauge theory, gauge fields are represented by group ele- 
ments rather than by fields, and the action is a periodic function of a gauge- 
invariant lattice field strength. The periodicity of the action entails spurious 
vacua. The principal advantage of noncompact actions, in which gauge fields 
are represented by fields, is that they avoid multiple vacua. 

Palumbo, Polikarpov, and Veselov Jj]] carried out the first gauge-invariant 
noncompact simulations. They saw a confinement signal. Their action con- 
tains five terms, constructed from two invariants, and involves (noncompact) 
auxiliary fields and an adjustable parameter. 

The present paper describes a test of a new way [0 of performing gauge- 
invariant noncompact simulations. Our action, which is similar to one term of 
Palumbo's action, is exactly invariant under compact gauge transformations, 
is a natural discretization of the classical Yang-Mills action, and reduces to 
Wilson's action when the gauge fields are compactified. In our version of 
Palumbo's method, there are fewer auxiliary fields, and they are compact 
group elements related to gauge transformations. 

We have used this method to simulate the U(l), SU(2), and £77(3) gauge 
theories on 8 4 and 12 4 lattices. For U(l) our Creutz ratios of Wilson loops 
agree with the exact ratios of the free continuum theory after a renormaliza- 
tion of the charge when the inverse coupling ft exceeds 0.5. Our SU(2) and 
SU (3) Creutz ratios clearly show quark confinement at ft = 0.5 and ft = 2, 
respectively. At much weaker coupling, our SU(2) and SU{3) Creutz ratios 
approximate those of continuum perturbation theory when the coupling con- 
stant is suitably renormalized. For SU(3) there is a scaling in the transition 
region 2.2 < ft < 3 with a lattice QCD parameter A L « 130 ±18 MeV, which 
is to be compared with the continuum value « 210 MeV. 



2 The Method 

For massless fermions, the continuum action density is ipij^d^ip. A suitable 
discretization of this quantity is iip (71)7^(11 + e M ) — i/j(n)]/a in which n is 
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a four-vector of integers representing an arbitrary vertex of the lattice, is 
a unit vector in the fith direction, and a is the lattice spacing. The product 
of Fermi fields at the same point is gauge invariant as it stands. The other 
product of Fermi fields becomes gauge invariant if we insert a matrix A^n) 
of gauge fields 

$(n)ip[l + igaA ll (n)]'il)(7i + e li ) (1) 

that transforms appropriately. Under a gauge transformation represented by 
the group elements U(n) and U(n + e M ), the required response is 

1 + iagA'^n) = U(n)[l + iagA^n)]U-\n + e M ) (2) 

or equivalent ly 

A'^n) = U{n)A^n)U-\n + + —U(n) \U~\n) - U'^n + e M )l . (3) 

ag 

Under this gauge transformation, the lattice field strength 

F^(n) = -[A„(n + e v )-A„(n)]--[A v {n + e (t )-A v (n)] 

+ ig [A v (n)Ar{n + e v ) - A^(n)A v (n + e„)] , (4) 

which reduces to the continuum Yang-Mills field strength in the limit a — > 0, 
transforms as 

KM = U(n)F^(n)U-\n + e, + e v ). (5) 

The field strength F^n) is antisymmetric in the indices // and u, but it 
is not hermitian. To make a positive plaquette action density, we use the 
Hilbert- Schmidt norm of F^(n) 

S = ^Tr[Fl(n)F, u (n)}, (6) 

in which the generators T a of the gauge group are normalized as Tr(T a T{,) = 
k5 a b- Because F^in) transforms covariantly (|5|), this action density is exactly 
invariant under the noncompact gauge transformation (|^). 

The gauge transformation (0) with group element U(n) = exp(—iagu a T a ) 
typically maps the matrix of gauge fields A^in) = T a A a ^{n) outside the Lie 
algebra, apart from terms of lowest (zeroth) order in the lattice spacing a. 
We accept this larger space of matrices. We use the action (0) in which the 
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field strength (|j) is defined in terms of gauge-field matrices A^in) that are 
the images under arbitrary gauge transformations of matrices A ^ (n) of gauge 
fields defined in the usual way, A°(n) = T a A^°(n) where the fields A a ^°{n) 
are real: 

1 + lagA^n) = V,(n)[l + iagAl(n)\W; l {n + e M ). (7) 

These gauge transformations are applied separately and link-wise to the 
gauge-field matrices A^(n). In general the group elements V^(n) and W ll (n + 
e M ) associated with the gauge field A^{n) are unrelated to those associated 
with the neighboring gauge fields A^{n + e u ), A u (n), and A u {n + e M ). But 
when the matrices and are equal and independent of direction 

V„(n) = W^n) = U(n) (8) 

for every vertex n, then they constitute a gauge transformation. 

Actually there is only one independent group element associated with 
each link. For by writing the equation (^) in the form 

1 + iagA^n) = V„{n)[\ + lagA^n^V-^n) 

xV^W-^n + e,), (9) 

we see that the departure of the gauge- field matrix A^n) from the Lie algebra 
is entirely due to the last matrix product V r At (n)H / / 7 1 (n + e M ). Thus there 
are as many auxiliary fields in this method as there are generators of the 
gauge group. For SU(N), the product VW' 1 has iV 2 — 1 generators, so our 
method involves iV 2 — 1 auxiliary fields in this case. We have not tried to 
parameterize these products and have instead accepted whatever matrices 
the link-wise gauge transformation (0) generated. 

Palumbo, in his version of this method [I], proposed a parameterization 
for the gauge-field matrices A^ that is slightly more general than is necessary 
and that leads to two more auxiliary fields than in our version. For instance, 
in the case of SU(N) his procedure involves iV 2 + 1 auxiliary fields which 
support a U(N) gauge invariance. His action also has a special term that 
may be needed to suppress these extra auxiliary fields. 

The auxiliary fields that describe the group element VW~ l are the prin- 
cipal defect of the present method. Their presence may be detected by 
measuring the average value of the path-ordered products of the factors 
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V^{n)W l {n + e M ) around the plaquettes of the lattice. For when the gauge- 
transformation condition @ is satisfied, the product around each plaquette 

P^{n) = V^W^in + e^V^n + e^W-^n + e^ + e^) 

x [V^n + e v )W~\n + e, + e v )^[V v {n)W;\n + e„)]t (10) 

is unity. 

To estimate the effects of the auxiliary fields, we have measured the mean 
values of these products 

0Ll n,(w 

at various values of (3 in our SU(3) simulations. We found that the auxiliary 
fields do reduce to gauge transformations in the continuum limit, /3 — > oo, 
but slowly: at (3 = 2, (P) = 0.18; at /3 = 3, (P) = 0.30; at /3 = 100, 
(P) = 0.71; and at = 1000, (P) = 0.87. 

The quantity 1 + igaA^n) is not an element L^in) of the gauge group, 
except when the real gauge fields A^°(n) all vanish. But if one compactified 
the fields by requiring 1 + igaA^n) to be an element of the gauge group, 
then the matrix A^ri) of gauge fields would be related to the link L^n) by 
Afj,(n) = (Lfj,(n) — l)/(iga), and the action @ defined in terms of the field 
strength (f|) would be, mirabile dictu, Wilson's action: 

_ k - MTrL»L„(7i + e M )L]> + e„)4(n) 
6 " Wfc • {U) 



3 Results 



We have tested this method by applying it to the U(l), ££7(2), and SU(3) 
gauge theories on 8 4 and 12 4 lattices. We used a multi-hit Metropolis algo- 
rithm with a flat measure for the real gauge fields A^°(n) and the relevant 
Haar measure for the group elements V^(n) and W^n). In most of our ini- 
tial configurations, the unitary matrices V and W and the hermitian gauge 
fields A ^ were randomized. For thermalization we allowed 50,000 sweeps 
for £7(1), 10,000 for SU(2), and 100,000 for SU(3). We allowed at least 20 
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sweeps between measurements. Our Wilson loops are ensemble averages of 
ordered products of the binomials 1 + iagA^in) rather than of the exponen- 
tials expliagA^n)} around the loop. 

For i7(l) and for /3 = 1/g 2 > 0.5, our measured Creutz ratios [Q] of Wilson 
loops agree with the exact ones of the free continuum theory apart from finite- 
size effects and a renormalization of the charge. For instance at f3 — 0.5 on 
the 12 4 lattice, we found x(2,2) = 0.147(1), x(2, 3) = 0.103(1), x(2,4) = 
0.090(1), x(3,3) = 0.049(1), x(3,4) = 0.034(1), and x(4,4) = 0.020(2). At 
the value (3 r = 0.45 of the renormalized charge, the exact Creutz ratios of 
the continuum theory are: x(2, 2) = 0.147, x(2, 3) = 0.102, x(2,4) = 0.093, 
X(3,3) = 0.048, x(3,4) = 0.037, and x(4,4) = 0.024. The agreement is 
excellent apart from the two largest loops which are too small by 8% and by 
17%, respectively, largely due to finite-size effects. 

At stronger coupling, the extra terms 

ig [A v (n)A^n + e u ) - A^n)A v (n + e„)] (13) 

in the lattice field strength F^in) eventually do produce a confinement sig- 
nal. For example, at (3 — 0.375, our measured Creutz ratios on the 12 4 
lattice are: x(2, 2) = 0.906(5), x(2, 3) = 0.909(21), x(2,4) = 0.85(10), 
X (3, 3) = 0.62(24), and X (3, 4) = 0.6(16). 

For SU{2) on the 8 4 lattice at /3 = 4/g 2 = 0.5, we found x(2, 2) = 
0.835(3), x(2,3) = 0.852(12), x (2,4) = 0.865(60), and x(3,3) = 0.94(23) 
which within the limited statistics clearly exhibit confinement. 

To test this method at weaker coupling, we compared our Creutz ratios 
with those of continuum perturbation theory. The tree-level perturbative 
formula for SU(n) may be expressed in terms of the function u(i,j) 

i i j j i ( 1 1 \ /., ,\ 

u \ l -,j) — ~ arctan — I — arctan log — H — - (14) 

3 i \i 2 j 2 J 

as 

n 2 - 1 

Xihj) = 27r2 p - u(i - 1, j - 1) + u(i,j - 1) + u(i - 1, j)] . (15) 

For SU(2) at (3 = 1, our six Creutz ratios x(2,2) = 0.1087(5), x(2,3) = 
0.0783(5), x(2,4) = 0.0696(5), x(3,3) = 0.0414(7), x(3,4) = 0.0299(8), and 
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x(4, 4) = 0.0189(10) are close to the x' s of the perturbative formula (|T5]) 
at the renormalized value of (3 r = 1.75, to wit: x(2, 2) = 0.1133, x(2, 3) = 
0.0786, x(2,4) = 0.0721, x(3, 3) = 0.0374, x(3,4) = 0.0286, and x(4,4) = 
0.0187. 

For SU (3) at (3 = 6/g 2 = 2 on the 12 4 lattice, we found in ten independent 
runs x(2,2) = 0.838(1), x(2, 3) = 0.826(3), x(2,4) = 0.828(13), x(3, 3) = 
0.793(42), x(3,4) = 0.47(25), and x(4,4) = 1.2(86). Within the statistics, 
these results robustly exhibit confinement. 

At much weaker coupling, our ratios approximate those of the tree-level 
formula (O) of continuum perturbation theory apart from finite-size effects 
and after a renormalization of the coupling constant. At (3 = 4 on the 8 4 
lattice, for instance, we found in one run x(2, 2) = 0.0878(1), x(2, 3) = 
0.0622(1), x(2,4) = 0.0554(2), x (3, 3) = 0.0319(2), x (3,4) = 0.0228(6), and 
x(4, 4) = 0.0132(1); whereas the tree-level formula ([jlj) at the renormal- 
ized value of /3 r = 6.03 gives x(2,2) = 0.0878, x(2,3) = 0.0609, x(2,4) = 
0.0559, x(3,3) = 0.0290, x(M) = 0.0222, and x(4,4) = 0.0145. Sim- 
ilarly at ft = 100 we found in one run x(2,2) = 0.00477(1), x(2,3) = 
0.00325(1), x(2,4) = 0.00289(1), x(3,3) = 0.00150(1), x(3,4) = 0.00104(2), 
and x(4,4) = 0.00053(3); whereas the perturbative formula ( ]I5]) at (3 r = 111 
gives X (2,2) = 0.00476, X (2,3) = 0.00331, X (2,4) = 0.00303, X (3,3) = 
0.00157, x(3,4) = 0.00120, and x(4,4) = 0.00078. The better agreement for 
the smaller loops is a finite-size effect. 



4 Scaling 



We used an 8 4 lattice to study the scaling of the lattice spacing a with the 
coupling constant g for SU(3). The two- loop result for the dependence of 
the string tension era 2 upon the inverse coupling P is 



oa 



a 
A 



2 exp 



lir 2 P 102 , 

1 log I 

33 121 6 1 33 



7T 2 /T 



(16) 



We expect this scaling formula to hold in a transition region where per- 
turbation theory is still valid and where the quark- ant iquark static potential 
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V(r) is a linear combination of a confining linear potential and a Coulomb 
potential. At one end of this region, V(r) is mostly linear; at the other end 
it, is mostly Coulomb. We found a transition region between (3 = 2.2 and 
(3 = 3. In this region, we used the interpolation 

X(i,j) = I [(4 - (3)oa 2 + {(3- l)xo(i,j)] (17) 

in which the string tension aa 2 is given by the scaling formula flIS|) and the 
perturbative Creutz ratio Xo{hj) is given by the tree-level formula (|i~5|). 

In the figure we have plotted the Creutz ratios x(hj) that we measured 
on an 8 4 lattice for 1.5 < (3 < 3. We also have plotted the interpolative 
formula flT7|) for 2.2 < (3 < 3 for various values of a/A 2 from 8 to 14 as solid 
curves. Our 2x2 Creutz ratios fit the interpolation flTTj ) in the transition 
region 2.2 < (3 < 3 with a / A\ = 14. Our ratios *(2, 3) and x(2, 4) fit it 
with o j A\ = 12 and 11, respectively. Our ratios x(3, 3) fit it with a / ' A 2 L = 8. 
Altogether our x's fit the interpolation (|P7j) for 2.2 < (3 < 3 with cx/A 2 = 
11 ± 3. A string tension y/a w 420 MeV implies that A L w 130 ± 18 MeV, 
which is to be compared with the continuum value of A^ ~ 210 MeV and 
with the lattice parameter Alw ~ 7.9 MeV of Wilson's method. 

At stronger coupling where the transition to a linear potential is complete, 
our j)'s fit the scaling formula ([16]) without any Coulomb term for 1.9 < 
(3 < 2.1 if we set er/A| ~ 25.0 ± 4. The corresponding value of A^ is 
A L w 85 ± 7 MeV. 
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Figure Caption 

The SU(3) Creutz ratios xihj) an d the scaling predictions for the string 
tension era 2 are plotted against (3. The solid lines for 1.9 < f3 < 2.1 represent 
formula flT6|) for era 2 with cr/Af = 25 ± 4. The solid curves for 2.2 < (3 < 3.0 
are the interpolation flT7| ) between the tree-level formula (|T5]) for the Creutz 
ratios xihj) an d the scaling prediction (|16|) for era 2 with a/A| = 11 ± 3. 
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